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ABSTRACT 

The purpose of present paper is to find applications of Fractional Calculus approach in Fluid 
Mechanics. In this paper by generalizing the instability phenomenon in fluid flow through porous 
media with mean capillary pressure with transforming the problem into Fractional partial differential 
equation and solving it by applying Fractional Calculus and special functions. 

Keywords: Fluid flow through porous media, Laplace transform, Fourier sine transform, Mittag - 
Leffler function, Fox- Wright function, Fractional time derivative. 

1. INTRODUCTION AND MATHEMATICAL PREREQUISITES 

Using ideas of ordinary calculus, we can differentiate a function, say, f(x)=x to the I st or II nd 
order. We can also establish a meaning or some potential applications of the results. However, can 
we differentiate the same function to, say, the 1/2 order? Better still, can we establish a meaning or 
some potential applications of the results? We may not achieve that through ordinary calculus, but 
we may through fractional calculus — a more generalized form of calculus. Fractional Calculus is the 
branch of calculus that generalizes the derivative of a function to non-integer order. In other words 
the fractional calculus operators deal with integrals and derivatives of arbitrary (i.e. real or complex) 
order. The name "fractional calculus" is actually a misnomer; the designation, "integration and 
differentiation of arbitrary order" is more appropriate. The subject calculus independently discovered 
in 17 th century by Isaac Newton and Gottfried Wilhelm Leibnitz, the question raised by Leibnitz for 
a fractional derivative was an ongoing topic for more than three hundred years. For a long time, 
fractional calculus has been regarded as a pure mathematical realm without real applications. But, in 
recent decades, such a state of affairs has been changed. It has been found that fractional calculus can 
be useful and even powerful, and an outline of the simple history about fractional calculus.. The 
various researcher investigated their investigations dealing with the theory and applications of 
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fractional calculus Free shear flows are inhomogeneous flows with mean velocity gradients that 
develop in the absence of cigarette, and the buoyant jet issuing from an erupting volcano - all 
illustrate both the omnipresence of free turbulent shear flows and the range of scales of such flows in 
the natural environment. Examples of the multitude of engineering free shear flows are the wakes 
behind moving bodies and the exhausts from jet engines. Most combustion processes and many 
mixing processes involve turbulent free shear flows. Free shear flows in the real world are most often 
turbulent. The tendency of free shear flows to become and remain turbulent can be greatly modified 
by the presence of density gradients in the flow, especially if gravitational effects are also important. 
Free share flows deals with incompressible constant-density flows away from walls, which include 
shear layers, jets and wakes behind bodies. Hydrodynamic stability is of fundamental importance in 
fluid dynamics and is a well-established subject of scientific investigation that continues to draw 
great curiosity of the fluid mechanics community. Hydrodynamic instabilities of prototypical 
character are, for example, the Rayleigh-Benard, the Taylor-Couette, the Benard-Marangoni, the 
Rayleigh-Taylor, and the Kelvin-Helmholtz instabilities. Modeling of various instability mechanisms 
in biological and biomedical systems is currently a very active and rapidly developing area of 
research with important bio technological and medical applications (biofilm engineering, wound 
healing, etc.). The understanding of breaking symmetry in hemodynamics could have important 
consequences for vascular biology and diseases and its implication for vascular interventions 
(grafting, stenting, etc.). When in a porous medium filled with one fluid and another fluid is injected 
which is immiscible in nature in ordinary condition, then instability occurs in the flow depending 
upon viscosity difference in two flowing phases. When a fluid flow through porous medium 
displaced by another fluid of lesser viscosity then instead of regular displacement of whole front 
protuberance take place which shoot through the porous medium at a relatively high speed. This 
phenomenon is called fingering phenomenon (or instability phenomenon). Many researchers have 
studied this phenomenon with different point of view. Fractional calculus is now considered as a 
practical technique in many branches of science including physics (Oldham and Spainier [13]). A 
growing number of works in science and engineering deal with dynamical system described by 
fractional order equations that involve derivatives and integrals of non-integer order (Bensonet al. 
[2], Metzler and Klafter [9], Zaslavsky [24]). These new models are more adequate than the 
previously used Integer order models, because fractional order derivatives and integrals describe the 
memory and hereditary properties of different substances (Poddulony [14]). This is the most 
significant advantage of the fractional order models in comparison with integer order models, in 
which such effects are neglected. In the context of flow in porous media, fractional space derivatives 
model large motions through highly conductive layers or fractures, while fractional time derivatives 
describe particles that remain motionless for extended period of time (Meerscheart et al. [8]). 

The phenomenon of instability in polyphasic flow is playing very important role in the study 
of fluid flow through porous media in two ways viz. with capillary pressure and without capillary 
pressure. The statistical view point was studied by Scheidegger and Johnson [21], Bhathawala and 
Shama Parveen [3] considering instability phenomenon in porous media without mean Capillary 
pressure. Verma [23] has also studied the behavior of instability in a displacement process through 
heterogeneous porous media and existence and uniqueness of solution of the problem was discussed 
by Atkinson and Peletier [1]. El-Shahed and Salem [11, 12] have used the fractional calculus 
approach in fluid dynamics, which has been described by fractional partial differential equation and 
the exact solution of these equations, have been obtained by using the Laplace transform, Fourier 
transform. Flow in a porous medium is described by Darcy's Law (El-Shahed and Salem [11]) which 
relates the movement of fluid to the pressure gradients acting on a parcel of fluid. Darcy's Law is 
based on a series of experiments by Henry Darcy in the mid- 19th century showing that the flow 
through a porous medium is linearly proportional to the applied pressure gradient and inversely 
proportional to the viscosity of the fluid. In one dimension, q represents "mass flow rate by unit area" 
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and is defined as, 

_ K dP 
Q ~ S dx 

where K is permeability, a parameter intrinsic to the porous network. The unit of permeability K is — 

9 1 9 1 dP 

. 8 is the kinematics viscosity has dimension L T , e.g., cm sec and — is the non-hydrostatic 
part of pressure gradient has dimension M LT 2 T~ 2 e.g., g cmT 2 sec~ 2 . thus the mass flow rate by 
unite area(q) has dimension c ™ - — $ — - = g „ Here, we considered homogeneous dimensions 
but in fractional calculus dimensions are inhomogeneous. 

Mathematical Prerequisites 

The Laplace Transform (Sneddon [20]) is defined as, 
L {/(*)} 

oo 

= j e~ st f(t)dt {Re (.s) 
o 

> 0) (1.1) 

Fourier transform of Weyl fractional derivative D y x f{x. t) is given by (Metzler and 
Klafter,2000,p 59A 12,2004) 

F{_ m Dlf(x,t)}= \k\r<p(k,t) 

The Fourier sine transform (Debnath [5]) is given as, 

u (n, t) 

oo 

2 

0 

The Error Function ([16]) of x is given as 
er f(x) 

X 

2 



uu 

2 f 

= — u(x, t) sinnx dx (1-2) 

VTT J 



= -j exp(-t 2 )dt (1.3) 



o 

And the complimentary error function of x is given as 
er f c {x) = 

2 -f"ex V (-t 2 )dt (1.4) 

In 1903, the Swedish mathematician Gosta Mittag-Leffler [10] introduced the function E a (z) which 
is given as — 

oo 

E » (z) = Ir(JTi)- W 

71 = 0 

Where z is a complex variable and T(a) is a gamma function of a. The Mittag-Leffler function is 

direct generalization of the exponential function to which it reduces for a = 1. For 0 < a < 1, E a (z) 

i 

interpolates between the pure exponential and a hyper geometric function — — . Its significance is 
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realized during the last two 100 years due to its involvement in the problems of physics, chemistry, 
biology, engineering and applied sciences. Mittag-Leffler function naturally occurs as the solution of 
fractional order differential equation or fractional order integral equations. The generalization of 
E a (z) was studied by Wiman [23] in 1905 and defined the function as 

oo 

E«,p(z) =X r(an + /?) ' (a ' p G C; Re(a) > 0 ' Re( ^ 

71 = 0 

> 0) (1.6) 

Which is known as Wiman' s function or generalized Mittag-Leffler function as Ea, x (z )=E„ (z ). 
The Laplace transform of (1.6) takes in the following form (Shukla and Prajapati [20]) 



e st t aj+f3-i £ W( xt «) dt 



j ,,-sf j-aj + P-1 £.(/> 
0 

(1.7) 



(s a - x)J +1 
Where E^{Z) = £j E ufi (z). 

The Fox-Wright function (Craven and Csordas [4]) is given as under 

V oo n^rcajk + bj) ^ 
- 2^=o n P =ir(c . k + dk) k[ , 

Where T(x) denotes the Gamma function of x and p and q are nonnegative integers. If we set bj = 1 (j 
= 1, 2, 3, p) 

and dj = 1 (j = 1, 2, 3, q ) then (1.8) reduces to the familiar generalized hyper geometric function 
(Craven and Csordas [4]). 

AC* <h» . = Z«--°Wtj$-S. <W 

For the study of generalized Navier - Stokes equations, El - Shahed and Salem [12] used the very 
special case of (1.1), given as 

oo 

j = 0 ' 

The Laplace Transforms of (1.10) is given by 



/ 



e~ st w{a,P; t) dt 



~s Ea 'P is) 



(1.11) 

The relationship between the Wright function and the Complementary Error function is given as, 



w 



(-i l;z) = er/ c (z/2) (1.12) 
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Riemann-Liouville fractional integrals of order u (Khan and Abukhammash [7]) is given as 
Let f (x) G L(a, b), u G C (Re(u) > 0) then 

</(x) = /£+/(*) 
1 f x f(t) 

= fool J^Wi dt (x > a) (113) 

is called R-L left- sided fractional integral of order U. 
Let f (x) G L (a, b), u G C (Re(u) > 0) then 

l r b f(t) 

- ml (t - dt (x 

< b) (1.14) 
is called R-L right-sided fractional integral of order u. 

The Laplace Transform of the fractional derivative (El - Shahed and Salem [12]) is defined by 

°° n-l 

j e~ st D a t f(t)dt= s a f(s) -^s"-'- 1 D a t f(PXn- 1< a < n) (1.15) 
0 ;'=o 

Theorem (Asymptotic expansion of Wiman function E a p(z )): Let 0 < a < 1 and P be an arbitrary 

complex number 

then 

/ I\ H 
^ exp I G« I G a 

2cr7ri J G — z 



E «^ z = ^f — "-^z — d G a-") 



We also use following integral (El - Shahed and Salem [12]) in terms of Wright function as, 

00 

J nsmnxE aia+1 (-n 2 Ct a )dn= ^—W -(^,1;-=) (1-17) 



2. STATEMENT OF THE PROBLEM 



If water is injected into oil saturated porous medium, then as a result perturbation (instability) 
occurs and develops the finger flow (Scheidegger and Johnson [19]). In this paper, our aim is to 
study one dimensional flow, x-indicating the direction of fluid flow with the origin at the surface, due 
to presence of large quantity of water at x = 0. We assume that water saturation at x = 0 is almost 
equal to one i.e.l and water saturation remain constant during the displacement process. Our picky 
interest in this paper is to look at the possibilities of transforming the problem in form of fractional 
partial differential equation with proper initial and boundary conditions. 



3. FORMATION OF THE PROBLEM 



The seepage velocity of water (Vw) and oil (8w) are given by Darcy's law [17] as 
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= -^K^ (3.1) 



5 W dx 



K 0 „ dP 0 



V 0 



And equation of continuity 

dS w dV w 

0 h 

dt dx 

= 0 (3.3) 

dS 0 dV 0 

0 h 

dt dx 

= 0 (3.4) 

Where K is the permeability of the homogeneous medium, Kw and Ko are the relative 
permeability of the water and oil, Sw and So are saturation of water and oil respectively, Pw and Po 
are the pressure in water and oil, phases 5w and 8o are the kinematics viscosities of water and oil 
respectively and (p is the porosity of medium. 

For inhomogeneous dimensions, considering the Time-fractional partial differential equations 
of the two phases as under: 

d a S w dV w 

0 1F :+ 17= °- (0<a<1) (3 ' 5) 
0 ^ + a7= 0 ' C0<«<U (3-6) 

For a = 1, equations (3.5) and (3.6) reduce to equations of continuity (3.3) and (3.4) respectively and 
from the definition of phase saturation [17], we have 

S w + S 0 = 1 (3.7) 

The capillary pressure PC is defined as pressure discontinuity between the flowing phases across 
their common interface and assumes the function of the phase saturation is a continuous linear 
functional relation as 

P c = PS W (3.8) 
Pc = Po ~ P w , (3.9) 

where [3 is constant. 

Relationship between phase saturation and relative permeability [18] is given by 

K w — S w 



K 0 = l-S w 3.10 

= s 0 

4. FORMATION OF FRACTIONAL PARTIAL DIFFERENTIAL EQUATION 

Now, we put the values of V w and V 0 (from (3.1) and (3.2)) in (3.5) and (3.6) respectively, we get 
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0 



d a S w = d {K w 
dt a ~ dx [S w 



dt 

d a S, 



K 



dt a dx {S 0 dx) 



dx 
dP r 



r) P 

Eliminating — ^ from (4.1) and (3.9), 

ox 



(0 < a < 1 ) 
(0 < a < 1 ) 



0 ^ = JM^(«)1 

dt a ox t 8 U/ \ ox ox I) 



(4.1) 
(4.2) 



(4.3) 



From (4.2), (4.3) and (3.7) we obtained 

A \ K f^ + ^]^_ ^L K d J± 
dx I \8 U , 8 n ) dx 8,,, dx. 



= 0 



(4.4) 



Integrating (4.4) with respect to x, we get 

K f*w + M^_ *w K dPc = _ B 
l8 w 8 0 ) dx 8 W dx 
Where B is the constant of integration, whose value can be determined. 
Equation (4.5) can be written as 



dP^ 
dx 



-B 



8 W I K w 8 0 ] 



+ 



d_Pc_ 
dx 



Ku, 8 n 



dP 0 



Substituting the value of — — from (4.6) in (4.3), we ar 



d a S d 
0_^ 0_ 

dt a dx 



K KodP L 
8 0 dx 



1 + 



o dx | B 

Ko_§w_ i , f^o §w) 

K w 8 0 \K W 8 0 y 



= 0 



(4.5) 



(4.6) 



(4.7) 



Pressure of oil (P 0 ) can be written as, 

Po =\(Po+P W )+ \{Po~P w ) 

- 1 

= p + -p c 

2 c 

where P is the mean pressure, which is constant. 
From (4.5) and (4.8) we get, 

K (K w K n ^ dP r 



B= - p5L_lzl °J± 
2 \S U , 8 n ) dx 



(4.8) 



(4.9) 



Substituting (4.9) in (4.7), we get 



0 



d a S„ 



+ 



ld_( dP r dS v 



dt a 2 dx I 8... dS,., dx 



0 (0 < a < 1 ) 



(4.10) 
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Taking K — — — = —A then (4. 10) reduces in the form, 

S w 3S W 



d 2 S w _ 1 d a S w 
dx 2 c dt a 



(4.11) 



Where C = — 

2P 

Equation (4.11) is the desired fractional partial differential equation of motion for water 
saturation, which governed by the flow of two immiscible phases in a homogenous porous medium 
and appropriate initial and boundary conditions are associated with the description as to 

S w (x,0) = 0, 
5 W (0,T) = 5 Wo 

< 1, (4.12) 
lim S w (x, T) — 0; 0 < x < co 

X— >oo 



5. FORMULATION AND SOLUTION OF PROBLEM 

Here we are formulating the generalized fractional partial differential equation by 
generalizing of Prajapati's et. al.[16] model equation which is given as it is in equation (4.1 1) 
From (4.11), we have 

d a S d 2 S 
dt a ~ dx 2 

Now, we generalizing the above equation by using fractional calculus operators in following form 



d a S w (x,t) a dPs w (x,t) 



dt a dtP 

Or 



o Z) t a S w (x, t) + a o D[S w (x, t) 

= C _oo Dl S w (x, t) + (p(x, t) (5.2) 

Where 0 Df is a Riemann-Liouville Fractional Derivative , ^ m D y x is a Weyl Fractional 
derivative^ (x, t)is a constant which describes the nonlinearity in the system. 

If we put a = 0, y = 2 and <p(x, t) = it converts into Prajapati's et. al.[16] paper 

d a S w (x,t) dPSJjcj) dVS w {x,t) 

ha „ = C +(p(x,t) (5.3) 

dt a dtP dxY rK J v J 

Solution-: Applying Laplace transformation both sides in equation (5.2), we get 
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n-1 ri-1 



s a S w (x,s) - ^ s B " r - 1 D t S w (x,o) + asPS w (x,s) - ^ as^' 1 D t S w (x,o) 

r=0 r=0 

= C_ m D Y t S w (x, s) + T(x7s) (5.4) 

Apply boundary conditions from (4.12) 

s a S w (x,s) + as? S w (x,s) = C_ooD[ S w (x,s) +f(x,s) (5.5) 

As is usual, It is convenient to employ the symbol - to indicate the Laplace transform with respect to 
the variable t. Also f(x, s) is Laplace transform of q)(x, t) [17] . 

Now we apply the Fourier transform with respect to the space variable x to the above equation both 
sides 

s a S^\K,s) + asPS^\K,s) 

= -C\K \VS W (K,s) + f(k,s) ( 5.6) 



(5 -ZL 

S w \K,s) 



5 "*^< S ) = ^Brt where b=<W 



f{k,s) 

(5.8) 



s a + asP +b 



Applying Inverse Laplace transform both sides 

s^K.t) = 2- 0 {-ay^<p\k,t-e)e a+ ^-^E r a %_ p) r+a (-bs a )ds 

Where b=C|/f| y (5.9) 

Now, Applying Inverse Fourier transforms both sides, we get the desired result, 

OO 

\r rt 



S w (x,t) = ^^T"J £ aHa-P)r-l j Q -ikx <p *( fcjt 
r=0 0 -00 



£ ) E a,a+(a-P) K-be a ) dfcrf£ ( 5 - 9 ) 

Which is a new and generalized result of Prajapati's et. al.[16] results. 

Now if we put a=0, y = 2 and cp — 0 

It converts into Prajapati's et. al.[16] results. 

This is easy to write in the form of wright functions 

5 w (x,t)= s w w(-^,l;-^) (5.10) 
On setting a = 1 and using (1.12), (5.7) reduces to, 

S w {x, t) = s Wo er f c {x I2<Q t) (5.1 1) 
Which is same as Prajapati's et. al.[16] results. 
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CONCLUSION 

The constitutive relationship model of generalization of the instability phenomenon in fluid 
flow through porous media with mean capillary pressure by applying fractional calculus is obtained. 
The exact solution of the generalized fractional partial differential equation in terms of Wright 
function by means of Laplace transform, Fourier transform with proper initial and boundary 
conditions has been found. If we set a = 1 then equation (5.10) reduces to (5.11), this method 
certainly useful than conventional method as the conventional method derived only for a = 1 
(equations (3.3) and (3.4)) whose solution given by equation (5.11). While this fractional calculus 
together with Fourier and Laplace transforms method presented in this paper also applicable for 0 < 
a < 1 whose solution given by equation (5.10). 
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